Auto-oscillation threshold, narrow spectral lines, and line jitter in spin-torque 
oscillators based on MgO magnetic tunnel junctions 
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We demonstrate spin torque induced auto-oscillation in MgO-based magnetic tunnel junctions. 
At the generation threshold, we observe a strong line narrowing down to 6 MHz at 300K and 
a dramatic increase in oscillator power, yielding spectrally pure oscillations free of flicker noise. 
Setting the synthetic antiferromagnet into autooscillation requires the same current polarity as the 
one needed to switch the free layer magnetization. The induced auto-oscillations are observed even at 
zero applied field, which is believed to be the acoustic mode of the synthetic antiferromagnet. While 
the phase coherence of the auto-oscillation is of the order of microseconds, the power autocorrelation 
time is of the order of milliseconds and can be strongly infiuenced by the free layer dynamics. 

PACS numbers: 72.25.Pn, 75.40.Gb, 75.60.-d, 75.75. +a, 85.75.-d 



I. INTRODUCTION 

A spin-polarized current traversing a magnetic multi- 
layer can impart spin angular momentum onto the mag- 
netization, giving rise to an interaction commonly re- 
ferred to as "spin-transfer torques" (STT]P^. Among 
the various magnetic systems exhibiting STT effects - 
domain walls^, single films'^, point contact J^, vortice^^, 
granular solid^ -, two model systems have been stud- 
ied more extensively: nanopillars patterned from metallic 
multilayers^ or from magnetic tunnel junctions (MTJ]P^. 
Despite the very different nature of electronic transport 
in the two systems (i.e. ohmic vs tunneling), it is pre- 
dicted that spin- torque effect in magnetic tunnel junc- 
tions should allow for large-amplitude spin-waves to be 
excited, well above the thermal occupation of the spin 
wave populations at equilibriumf^. This phenomenon of 
current-driven magnetization oscillations is of wide in- 
terest for two reasons. First, such spin-torque oscilla- 
tors reveal some aspects of magnetization dynamics in 
the presence of a strong frequency nonlinearity, which 
leads to nontrivial modifications to the excitations. Sec- 
ond, these systems have the potential for applications as 
compact microwave oscillators for which the frequency 
is readily tunable with applied magnetic fields and cur- 
rents.^ 

Using spin-transfer torques to excite magnetization 
oscillations in spin-valve nanopillars is now a well- 
established experimental technique. The signature of 
genui ne au to-oscillation is a clear current threshold be- 
I^q^Iq^MSIB of the microwave signal with three concomi- 
tant features: a narrowing of the spectral line, a drastic 
increase of the line power, and a single mode spectrum 
with no flicker (1//) noise. Despite intense research ac- 
tivity, auto-oscillations that exhibit all these signatures 
together are still not commonly observed in tunnel junc- 
tions. The debate remains open as to whether the physics 
of magnetic tunnel junction places some intrinsic restric- 



tion on the range in which auto-oscillation can be ob- 
served, or rather if the difficulty arises from material 
limitations, for instance the possibility of fabricating tun- 
nel junctions that can withstand the high current densi- 
ties required to excite auto-oscillations. In comparison 
to spin- valves, three physical differences are anticipated 
in tunnel junctions. First, the spin- transfer torque may 
include a large component with a field-like symmetr}^^. 
Second, there may be a nontrivial voltage bias depen- 
dence, which would not allow sufficiently large currents 
to be obtained. Third, the voltage drop across the tun- 
nel barrier may induce significant heating }^ A number of 
experiments have been conducted recently with the aim 
of resolving some of these issues j -*^^ ' -*^^ ' ^^ ' ^-*^ ' but as far as 
we can ascertain, none has provided any clear signature 
of the auto-oscillatory regime. 

In this article, we present an experimental study 
of magnetization auto-oscillations in MgO-based tunnel 
junctions, driven by spin-transfer torques. We work with 
tunnel junctions that have been submitted to large cur- 
rent densities in order to stabilize the MgO barrier prop- 
erties. Upon applying increasing currents, we observe a 
strong narrowing of the spectral line, down to 6 MHz at 
an oscillation frequency of 6.06 GHz, which is accompa- 
nied by a drastic power increase at the emission thresh- 
old, yielding spectrally pure oscillations free of flicker 
noise. In contrast to the results presented in Ref. |2TJ 
the frequency dispersion of the observed auto-oscilation 
mode is consistent with that of the acoustic mode of the 
synthetic antiferromagnet (SAF). Unexpectedly, the elec- 
trons transferring spins to the free layer can destabilize 
both this free layer and the synthetic antiferromagnet for 
the same current polarity. Finally, we present some ev- 
idence of inhomogenous broadening of the spectral line 
when the linewidth is of the order of a few MHz; we show 
that the measured linewidth and lineshape is a result of 
a very narrow line that drifts faster than the accessi- 
ble capturing times. A spectral line envelope is therefore 
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measured in most cases, giving bounds for the power cor- 
relation time and the oscillation phase coherence time. 



II. SAMPLES AND EXPERIMENTAL 
METHODS 

Our nanopillars are fabricated from stacks of composi- 
tion PtMn (20) / CoToFeso (2) / Ru (0.8) / CoeoFesoBso 
(2) / [Mg 0.9 / natural oxide] / Co6oFe2oB2o (3), where 
figures in parentheses denote film thicknesses given in 
nm. The films were provided by the SINGULUS com- 
pany and grown using their Timaris sputter deposition 
cluster. The devices are rounded rectangles with lateral 
dimensions of 70x140 nm^ (Fig. IB, inset). Devices 
with other dimensions (50x100 and 100x200 nm^) also 
showed similar auto-oscillations, with minor differences 
to be reported elsewhere. For a chosen virgin device, 
the resistances of the P (parallel) and AP (antiparallel) 
states are 136 and 154 Q (magneto-resistance is 13%). 
The resistance- area product is 1.4 ^7 /im^. A current of 1 
mA (10^ A/cm^) generates an Oersted field of the or- 
der of 1.4 mT at the edge of the nanopillar. In our 
convention, positive currents describe a flow of electrons 
from the synthetic antiferromagnet layer to the free layer, 
which experimentally favors an AP orientation (Fig. IB). 
Material fatigue occurs for bias voltages above 0.4 V, 
after which we observe a large decrease in the junction 
resistance, with larger currents required to attain thresh- 
old for the auto-oscillations. However, we did not observe 
any qualitative changes to the magnetic properties of the 
free layer after degradation. 

Following the methods described in Ref. [23l the 
nanopillars are inserted in series in coplanar waveguide 
structures using design rules that allow an electrical 
bandwidth better than 50 GHz.^^ The voltage noise spec- 
tra were recorded under dc currents with a resolution 
bandwidth of 2 MHz. For a reason that will become ap- 
parent in section IV, we have recorded our spectra by 
averaging multiple frequency sweeps, except Fig. 3C in 
which a single sweep curve is presented. 



III. QUASI-STATIC PROPERTIES 

The quasi-static properties of a representative virgin 
nanopillar are presented in Fig. 1. The resistance versus 
easy axis field loop (Fig. la) indicate uniform magne- 
tization configurations in the remanent states. In our 
convention, positive fields are defined to favor an an- 
tiparallel orientation of the free layer with respect the 
top layer (reference layer) of the synthetic antiferromag- 
net (SAF). The SAF layers are not perfectly compen- 
sated, which results in a loop shift slightly favoring the 
P orientation. The Stoner-Wohlfarth astroid (Fig. 2B) 
confirms the uniaxial anisotropy expected from the de- 
vice elongated shape. The astroides exhibit the break 
that is usually present in those size ranges, which is 
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FIG. 1: Quasi- static properties of the nanopillars in the virgin 
state, (a) Resistance versus easy axis field minor loop under 
an applied current of 0.1 mA. Inset: high field curve, (b) 
Voltage versus applied current curve taken in zero applied 
field. Inset: scanning electron micrograph of the device after 
the pillar etching. 
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FIG. 2: (Color online) Stoner-Wohlfarth astroids of a device 
after the aging procedure using a bias current of 8 mA. A 
dc current is applied during each astroid measurement. The 
applied currents are -5 mA, 0.1 mA and 5 mA for the left, 
central, and right panels, respectively. 



representative of magnetization switching from P to AP 
through an incoherent micromagnetic path due to the 
inhomegeneity of the dipole field of the uncompensated 
SAF. The spin-flop transition of the SAF occurs near 100 
to 150 mT depending on the device (Fig. la, inset), i.e. 
in a region sufficiently far from the fields that we will 
employ to study spin-torque induced auto-oscillations. 

In the virgin state, the voltage versus applied current 
loops (Fig. IB ) reveal switching between the two resis- 
tance states at critical currents around ±1.5 mA. These 
switching thresholds increase with sample aging, which 
is accelerated through the use of large applied currents 
yielding voltages above 0.4V. The switching currents sta- 
bilize typically near 5 mA, provided that the bias never 
exceeds 8 mA. Expect for Figure 1, we have systemati- 
cally used this aging procedure to stabilize the junction 
properties. We emphasize already that these free layer 
switching thresholds are very different from the auto- 
oscillation threshold. 

The consequence of spin torque can also be observed 
by looking at the resulting distortion of the Stoner- 
Wohlfarth astroid of the free layer (Fig. 2). Under posi- 
tive currents, the spin-torque favors the AP configuration 
of the free and reference layers. In addition to a global 
shrinking which results from current-induced heating, the 
astroid becomes progressively distorted and a re-entrant 
border expands near the right corner of the astroid (Fig. 
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FIG. 3: Illustration of the current threshold behavior: Volt- 
age noise power spectra (linear scale, arbitrary units) versus 
current for a junction in the P state at an easy axis field of 
-4 mT. The device has been submitted to an aging procedure 
using a bias current of 8 mA. 



2C). This type of distortion is expected if the spin-torque 
has a pure Slonczewski form (see Fig. 9 in Ref. 26). Un- 
der negative currents (Fig. 2a), the current reduces the 
astroid in all directions by pushing the boundaries toward 
zero field, recalling that the current significantly heats 
the junction. In negative current, the astroid distortion 
is difficult to observe, indicating a lower spin torque per 
unit current. 



IV. SPIN TORQUE INDUCED 
AUTO-OSCILLATIONS 



In this section, we identify the current and field re- 
gions in which the applied spin-torques effectively favor 
auto-oscillation of some magnetic degrees of freedom of 
the nanopillars. We have analysed the voltage noise un- 
der currents up to 8 mA, with easy and hard axis fields 
up to 0.15 T, including above the spin- flop transition of 
the SAF. Thermal spin- waves are observed in most re- 
gions, but clear auto-oscillation is only observed at low 
field when the free layer and SAF layers are in the P 
configuration. In contrast to other studies,^^^the auto- 
oscillations we measure are also observed for zero applied 
external fields. Indeed, we believe this to be one of the 
first experimental demonstrations of field-free current- 
driven excitations in a magnetic tunnel junction for which 
a clear auto-oscillation threshold is observed. 

Let us now detail how the crossing of the threshold 
happens and how the auto-oscillations build up (Fig. 3). 
When the current is increased, a narrow line (F or 2F) 
suddenly emerges from the noise floor slightly below 3 
mA. The fundamental frequency (F) of this oscillation 
at zero field is 6 to 8 GHz, which varies from sample to 
sample and depends on sample aging. The frequency F 
is almost independent of the current (Fig. 3) until the 
oscillation ceases when the ground state of the nanopillar 
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FIG. 4: (Color online) Voltage noise signature of the au- 
tooscillation mode in various experimental conditions for a 
device after the aging procedure using a bias current of 8 mA. 
(a) Spectrum recorded in averaging mode from 0.1 GHz to 20 
GHz for easy axis field of -30 mT and 8 mA in the P state, (b) 
Narrow-band spectrum in averaging mode recorded for zero 
field and 3.4 mA in the P state, (c) Narrow-band spectrum 
of the second harmonic (2F) recorded in single sweep mode 
for an easy axis field of -23 mT and 4 mA in the P state. The 
quality factor for 2F is Q = 2600. 



switches to AP slightly below 8 mA (not shown). 

In Fig. 4, we present the spectral signature of auto- 
oscillations in several representative experimental con- 
ditions. In Fig. 4a, we present the broadband spectrum 
across the frequency range of 0.1-20 GHz which shows the 
absence of any significant 1// noise. Apart from higher- 
order harmonics, there is no other feature in the whole 
0.1 to 20 GHz interval, indicating that the ground state 
and the excited mode are uniquely defined. In Fig. 4b, 
a narrower band measurement is presented for an auto- 
oscillation at exactly zero applied field and under a cur- 
rent right at the emission threshold. Finally in Fig. 4c, 
we present the narrowest spectral line we have observed 
at room temperature, which was obtained after optimiz- 
ing the field magnitude and its orientation in the film 
plane. The spectral line was obtained using a single-shot 
measurement to avoid any envelope effects due to fre- 
quency jitter and inhomogeneous line broadening. Note 
that it is the second harmonic (2F) that is presented in 
Fig. 4c, as no signal at F was observed above the auto- 
oscillation threshold at that applied field. 

In the auto-oscillation regime, we observed noticeable 
jitter in the spectral line position. At constant applied 
currents and fields the measured line shape can vary be- 
tween each measurement, where the linewidth can fluc- 
tuate, for example, from 3 MHz as shown in Fig. 4c, to 
18 MHz in a worse case. This fluctuation suggests there 
is a significant jitter in the oscillation frequency during 
the acquisition of the spectral data. To get data with 
statistical significance, the spectra must be recorded in 
averaging mode, in order to ensure that the line fluctu- 
ation translates into a stationary line envelope. Under 
such measurement constraints, the best measured line 
envelope width was not better than 6 MHz, correspond- 
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FIG. 5: (Color online) Frequency and width of spectral line 
envelope, as a function of applied current at zero applied field. 
The device has been submitted to an aging procedure using 
a bias current of 8 mA. 



ing to a quality factor of 1000. Above threshold, the 
width of the line envelope decreases to a flat minimum in 
general (Fig. 5), until it increases with line shapes start- 
ing to more significantly vary from one spectrum to the 
next, indicating either unstable ground state or unstable 
oscillation. 

A prerequisite for observing the auto-oscillations in ex- 
periment is an adequate preparation of the initial ground 
state. In Fig. 6, we present the power spectra in a density 
map as a function of applied easy axis field for an applied 
current of 4 and 6 mA, which are above threshold. Both 
increasing and decreasing branches of the field sweep are 
presented. The auto-oscillation is observed only in the P 
state. While the second harmonic (2F) is always present, 
the fundamental (F) and third harmonic (3F) are only 
seen in a restricted field interval. We note that the pres- 
ence of the fundamental F is always observed systemat- 
ically with the third harmonic 3F (Fig. 6B). We also 
note that F and 3F are observed only in the field region 
where the free layer magnetization is not fully saturated 
(see Fig. lA). This indicates that the appearance of F 
and 3F is not related to a change of the oscillation mode, 
but rather corresponds to a change of the sensitity of 
the magneto-resistance response to the (unaltered) oscil- 
lation mode. The width of the line envelope appears to 
be independent of the easy-axis field in the monostable 
P area, i.e. with typical variations between 8 to 19 MHz 
for fields in the range of -40 to -20 mT (Fig. 6A). At 
high negative fields, this linewidth floor comes with the 
absence of the power at F and 3F (Fig. 6B). 

In addition to the auto-oscillation, a faint lower fre- 
quency mode with a large linewidth (400 MHz) can gen- 
erally be detected just before the AP to P (Fig. 6C, 
mode labelled "free layer") or the P to AP (not shown) 
transition. An extrapolation of the frequency of this faint 
mode indicates that it approaches zero at the anisotropy 
field of the free layer. This suggests that this mode cor- 
responds to the quasi-uniform spin wave mode of the free 
layer. The current dependence of the power of this peak 
(not shown) confirms that this mode is thermally popu- 
lated rather than being a persistant, spin-torque induced 
auto-oscillation. 

To explore further the nature of the excited mode at 
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FIG. 6: Voltage noise properties as a function of the easy axis 
field and ground state. The device has been submitted to an 
aging procedure using a bias current of 8 mA. In panels A 
and B the current is 6 mA and the field is first set negative 
to prepare the P state. Panel A: line envelope width of the 
2F peak. Panel B: noise power densities (log scale). In panel 
C, the current is 4 mA and the field is first set positive to 
prepare the AP state. 
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FIG. 7: Voltage noise spectra (log scale, arbitrary units) as a 
function of increasing hard axis field for a junction biased with 
+6 mA and -4 mT along the easy axis at 300K for a junction 
prepared in the P state. The device has been submitted to 
an aging procedure using a bias current of 8 mA. 



zero field, we varied the hard axis field at constant ap- 
plied current and constant easy axis field (Fig. 7). When 
far from the center of the Stoner-Wohlfarth astroid, the 
width of the line envelope is observed to be very broad 
(300 MHz), which leads us to speculate whether the ob- 
served line corresponds to one mode or several modes 
with closely-spaced frequencies (see Fig. 8C for hard axis 
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FIG. 8: Variation of the (a) spectral line power, (b) width of 
the line envelope, and the (c) frequency, which are extracted 
from Lorentzian fits to the power spectrum for a junction 
biased with +6 mA and -4 mT along the easy axis at 300K. 
The device has been submitted to an aging procedure using 
a bias current of 8 mA. 

fields negative and greater than -2 mT). The width of the 
line envelope reduces to 6 MHz as the hard axis field ap- 
proaches the middle axis of the astroid (Fig. SB). Note 
that this minimum in the width of the line envelope cor- 
relates with a zero slope of the center frequency at the 
same applied hard axis field. 

This drastic decrease in the width of the line envelope, 
which in some cases represents a reduction by a factor 
of 50, is accompanied by a large increase of the peak 
in the noise power density in such a manner that the 
total power of the auto-oscillation is almost independent 
of the hard axis applied field (Fig. 8A). This is another 
indication that the spectral line we measure corresponds 
to the envelope of a narrow line, whose power remains 
constant but suffers an excursion in the frequency within 
an interval that depends critically on the field. This is 
in stark contrast to the behavior observed in spin valves 
in which free layer excitations are measured^^. There, 
an increase of the hard axis field leads to a reduction in 
the linewidth due to an accompanying reduction in the 
frequency nonlinearity. 

V. DISCUSSION 

The auto-oscillation mode we observe at zero applied 
fields cannot correspond to an eigenmode of the free layer 
for the following reason. The frequency of a uniform pre- 
cession mode of the free layer decreases to zero, becom- 
ing soft at an applied field along the hard axis that is 
equal to the shape anisotropy field. From Fig. 7, it is 
clear that this is not the case: the mode frequency is al- 
most unchanged when the field approaches the hard axis 
corners of the astroid at -10 mT. A complementary ar- 
gument to demonstrate that our auto-oscillation is not 
a free layer mode can be made using the current depen- 
dence of the emission frequency. This dependence ex- 
hibits a very small blue shift (Fig. 5), while for a free 



layer mode there would be a strong red shift for in-plane 
applied fields, as expected from spin- wave theory. 

The frequency and its concave dependence on the hard 
axis field (Fig. 8C) recalls that of the acoustic branch of 
the SAF eigenexcitations. Reasonable numerical agree- 
ment is achieved (not shown) when assuming that the 
SAF layers have magnetizations of 0.9 T, that they are 
coupled by interlayer exchange coupling of -2.4x10"^ 
J /m^ , and that the exchange pinning strength is 90 mT 
Those numbers are based on our experience and 
are consistent with the experimental spin- flop fields. The 
interlayer and exchange pinning coupling strengths vary 
from sample to sample, and decrease with sample aging, 
probably because of thermal diffusion of Mn or Ru atoms 
during aging. 

Our observations of auto-oscillations under positive 
currents (electrons flowing from the SAF to the free layer) 
are unexpected. Indeed, it is well-established for metal- 
lic spin- valves that for a dc current that destabilizes the 
free layer, it stabilizes in turn the reference layer, and 
vice-versa. This is not the case here: a positive cur- 
rent destabilizes both the free and the SAF layers, by 
leading to a well-defined auto-oscillation of the synthetic 
antiferromagnet, but also switches the magnetization of 
the free layer at a slightly higher current (Figs. 1 and 
2). This behavior is suggestive of a non trivial bias 
dependence of the spin torque in these low resistance- 
area product magnetic tunnel junctions. One may won- 
der whether the spin torque solely originates from elec- 
trons crossing the oxide (resulting in free layer switching) 
or if it also partly originates from the diffusive trans- 
port in the synthetic antiferromagnet (resulting in SAF 
auto-oscillations). More elaborate transport calculations, 
based on realistic stacks in which the magnetization dy- 
namics of the free and synthetic antiferromagnetic layers 
are taken into account self-consistently, may be required 
to reproduce these experimental observations. Further- 
more, one would also need to account for realistic tunnel 
barriers, which are more "permeable" for electron trans- 
port than perfect MgO barriers. 

Another important issue concerns the origin of the 
line fluctuations, keeping in mind that the frequency 
and oscillator power are intimately related through the 
stro ng frequen cy nonlinearity with the oscillation ampli- 
tude.'SESESEUpj.Qj^ our findings, we shall distinguish two 
limiting cases: when the free layer is fixed, or when it un- 
dergoes substantial fluctuations. 

The first case corresponds to applied fields at the cen- 
ter of the Stoner-Wohlfarth astroid (Fig. 7), or to fields 
that further stabilize the free layer in a P ground state 
(Fig. 4C). In those cases, we are usually near a minimum 
of the width of the line envelope. We have observed fre- 
quency fluctuations that were slow enough to be visible 
from one spectrum measurement to next, while they were 
fast enough to be seen as random, and not as a gradual 
drift. This indicates that when the power and frequency 
deviate from their average, the excess (or default) power 
relaxes with a characteristic time scale that is compara- 
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ble to the time needed to record a line in a spectrum, 
i.e. of the order of some mihiseconds with our experi- 
mental set-up. In contrast, the phase coherence time of 
the oscillator is related to the best observed linewidths, 
and is therefore of the order of a few microseconds. This 
situation corresponds to the high temperature limit of 
Ref. ,32, in which the measured linewidth appears as an 
inhomogeneous broadening: power fluctuations are slow 
enough to result in a modulation of the oscillation fre- 
quency. This is corroborated by the fact that when the 
nonlinear ity vanishes (i.e. when experiment aly we have 
dF/dl = 0), the frequency is then expected to be immune 
to power fluctuations and the experimental line envelope 
width reaches indeed a minimum. 

In addition, we find that the width of the line envelope 
increases as the P-to-AP switching field is approached 
(Fig. 5), or when the free layer magnetization is tilted 
along the hard axis under a magnetic field (Figs. 7 and 
8). This behavior suggests that the free layer dynamics 
plays an important role on the auto-osc illation o f the SAF 
sub-system. However, present theories I ^^ I ^^ I ^Q I model the 
linewidth from noise sources that are solely internal to 
the oscillator, and consider this noise to be white. In 
our case, we believe the free layer dynamics is a signif- 
icant noise source for the SAF oscillations, either as a 
random field through the dipolar coupling, which would 
appear as a fluctuation in the oscillation frequency, or 
as a random modulation of the applied current through 
a varying spin-transfer efficiency, which would appear as 
a fluctuation of the supercriticality (I/Ic). In particu- 
lar, the free layer noise is peaked at the frequency of the 
lowest energy magnon (the main quasi-uniform ferromag- 
netic resonance mode). Since this mode goes soft when 
the switching field is approached, its population diverges, 
thereby contributing a strong colored noise to the torques 
driving the oscillations in the SAF. This may be the ori- 
gin of the increase in the observed linewidth near the free 
layer switching transition. 

VI. CONCLUSION 

In conclusion, we have identified the regions of the 
field-current parameter space for spin-torque induced 



auto-oscillations of the synthetic antiferromagnet layers 
in MgO-based tunnel junctions. The observed oscilla- 
tion mode exists in zero applied field and is the high 
amplitude counterpart of the acoustic eigenexcitation of 
the synthetic antiferromagnet layers. After sample ag- 
ing, the prerequisites for auto-oscillation include a very 
low resistance- area product, a ground state in the paral- 
lel configuration, and a weak or vanishing hard axis field. 
Upon increasing the current, three signatures of an auto- 
oscillation threshold are observed: a drastic spectral line 
narrowing down to 6 MHz (for an oscillation frequency of 
6.06 GHz), a substantial increase in the oscillation power, 
and a spectrally pure spectrum that is free of 1// noise. 
In contrast to predictions from simple macrospin mod- 
els, the auto-oscillation of the synthetic antiferromagnet 
is observed for a current polarity that also destabilizes 
the free layer. This results suggests that more complex 
spin-transport phenomena need to be accounted for in 
simulations of magnetization dynamics in low resistance- 
area tunnel junctions made of realistic stacks. Finally, 
the emission frequency is observed to wiggle randomly 
during the measurement time, such that the measured 
line is the envelope of the line frequency excursion. We 
suggest that noise sources from other magnetic layers in 
the device, such as the free layer in our case, need to 
be accounted for in further quantitative comparisons be- 
tween theory and experiment. 
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